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Abstract 

Improving upon the previous treatment by Mueller and Son, we derive the Boltzmann equation 
that results from a classical scalar field theory. This is obtained by starting from the corresponding 
quantum field theory and taking the classical limit with particular emphasis on the path integral 
and perturbation theory. A previously overlooked Van-Vleck determinant is shown to control the 
tadpole type of self-energy that can still appear in the classical perturbation theory. Further 
comments on the validity of the approximations and possible applications are also given. 
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I. INTRODUCTION 



There are many situations in nature where classical field theory is the most direct and 
efficient way of describing the system under study. This happens when the occupation 
number becomes large enough that the quantum commutator becomes irrelevant. Such 
situations include most of macroscopic electro-magnetic systems and the dense systems 
created in colliding heavy ions. In studying heavy ion collisions, both the particle degrees 
of freedom and the classical field degrees of freedom becomes relevant: The final state of 
a collision consists mostly of the particle degrees of freedom whereas the initial state of 
iwo approaching nuclei can be efficiently described by a classical non-Abelian gauge field 
3, y| • A natural framework to encompass both elements is the coupled system of classical 
field equation and the Boltzmann type of kinetic equation with the mean fieldjj]. 

An important problem in this framework is the conversion of classical field degrees of 
freedom into particle degrees of freedom and their subsequent thermalization. In this respect, 
the equivalence of the classical field theory and the Boltzmann equation proposed by Mueller 
and Sonj^ is significant since it has a potential of providing a consistent framework for the 
thermalization in heavy ion collisions. The current paper improves on the work of Mueller 
and Son in the following aspects. 

First, going from quantum to classical many body theory almost always involves the 
Wigner function ^, 0, Q]- In Ref.yj, this point was overlooked. In this work, we show 
that the Wigner functional does appear in the formulation due to the fact that the density 
operator is in general non-diagonal. 

Second, the Feynman rules used in Ref . Q for the classical field are a mixture of quantum 
and classical ones: The vertex rule is classical but the propagators are quantum mechanical. 
We show that consistent formulation only involves classical vertices and propagators. 

Third, formulation of classical path integral involves a functional determinant related to 
the Jacobi field. Therefore, perturbation theory based on classical path integral involves 
ghost fermions Q. It turned out that the c.ass.ca. perturbat.ou theory as formulated 
in Ref.p still contains a quantum effect in the form of the tadpole self-energy. The ghost 
contribution is shown to compensate this remaining quantum effect, thus making the theory 
consistently classical. ^ 

Despite these gaps in the formulation, the final result in Ref. is essentially correct in 
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the / S> 1 limit. The goal of this paper is to provide a more consistent derivation of the 
kinetic equation corresponding to the classical scalar field starting from the quantum scalar 
field. We take particular care in treating the path integral and the perturbation theory in a 
consistent manner. We also comment on the validity of this classical method and possible 
use in the context of understanding the thermalization in heavy ion collisions. 



II. CLASSICAL PATH INTEGRAL FROM QUANTUM PATH INTEGRAL 

To be specific, consider a real scalar field theory defined by the following Hamiltonian 

The equations of motion are 

(2, 

(3) 

where we defined E{(j)) = V^</' — V'{(j)) for notational convenience. 

In many publications since 1980's E. Gozzi and his collaborators have been extensively 
studying the properties of classical path integrals (See Refs. and references therein.). 
The starting point in these studies is the following transition probabihty between two points 
in the phase space 

P[ff, Ti-J; tf\(fi, vrf ; ti] = 6[(ff - ^c{tf\^i, vrf )] ^[vrj - 7r^^(t/|¥?i, vrf )] (4) 

Here ^ and vr*^ are the generalized coordinate and momentum and and vr^ are the solutions 
of Eqs. ((21) and (jHI) with the boundary conditions y9c(ij|v^j, Trf ) = and 7r)f (tjlyjj, vrf ) = vrf. 
The evolution of the initial phase space density Pci[<^, tt*^; tj] is given by 

In this paper we use the square bracket notation J[d(f)] to indicate a functional integral at a 
fixed time and the curly D notation / Vcp = J Hili [d(l)i\ to indicate a functional integral over 
space and time. We also implicitly absorb any normalization constants into the definition 
of [d(f)]. 
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To turn Eq.ijSJ into a path integral, we use the fact that transition probabihties satisfy 

J dy P[z\y]P[y\x] = P[z\x] (6) 

Dividing the time between the final time t/ = ^at+i and the initial time ti = to into 
equal intervals we obtain 

. N N 
i=l j=0 

where we have suppressed the time arguments in P. 

For sufficiently small At = {tf —ti)/{N+l), we should be able to solve the classical equa- 
tions of motion approximately. The simplest finite difference method is the Euler method 
based on the following approximation of the time derivative 

~ ZM^Zi!^ + 0(At2) (Backward Euler) (8) 

or 

/(t,) ~ Zi^^±lLz/M + o(At2) (Forward Euler) (9) 
We use here the Backward Euler method for ip 

ipk+i = ipk + nt+iAt + OiAt^) (10) 
and the Forward Euler method for vr*^ 



TT 



fc+1 



7Tt + E{ipk)At + 0iAt^) (11) 



so that ipk+i can be expressed solely in terms of ipk and vr^. This is, of course, not the only 
choice. Different discretization method in general leads to different form of the discretized 
path integral although they are all equivalent in the At ^ limit. 
For small enough At then, 

Pb,+i, vrlVil^,, vrj] = - ^, - vr|V,At]5[vrJVi - vrj - E{if,)At] (12) 

With this form, it is easy to check that the preservation of probability 
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is trivially satisfied. Using dummy variables x ^i-nd vr^ to express the (5-functionals, we finally 
get the classical path integral for the evolution of the density functional 

N N+l I N 



j=0 k=l \ j=0 

. N N f ^ 

j IIKI n [dXk] exp -z E Xfc (vrr+i " vr^ - E{^,)M 



j=0 k=l \ j=0 

N N / N 

X 

1=0 k=0 \ k=0 

X Pel bo, 71";^; ^o] (14) 



How closely can we reproduce the classical path integral Eq. p4|) from the corresponding 
quantum path integral? To answer this question, we need to go back to the basics of 
deriving a quantum mechanical path integral. Given a Hamiltonian operator H, we would 
like to rewrite the matrix element of the evolution operator, (0j|f/(t/,tj)|0j), as a path 
integral. Dividing the time interval into many small segments and using the fact that 
U{t,t')U{t',t") = Uit,t"), we can write 

{(l)N+l\U(tN+l,to)\(t)o) = {(j)N+l\U(tN+l,tN)U(tN,tN-l) ■ ■ ■U{ti,to)\(t>o) (15) 

There are many ways of inserting the resolutions of identity to turn this expression into a 
path integral just as there are many ways of discretizing the classical equation of motion. 
The prescription which most closely resembles Eqs. ljlOj) and (jllll turned out to be inserting 

1 = J M0yt][(i7rfc]|0fc)(0fc|7rfc)(7rfc| (16) 

between U(tk+i,tk) and U(tk,tk-i) for all 1 < < and inserting 

1= J [rfvTjv+i] IvTAf+i) (tTjv+i | (17) 

between {(f)N+i\ and U{tN+i,tN)- In this way, we get 

{<pN+i\U{tN+i,to)\<po) = J lildM J IIM 



k=l k=l 

/ N+l N+l \ 

X exp i E TTk{(j)k - 0fc-i) - i H{iTk, (pk-ijh (18) 

\ k=l k=l / 

The time evolution of a density operator needs two such path integrals since the matrix 



element 



(0/|p(t/)|0/) = {<Pf\U{tf,U)p{U)U{U.tf)\^f) 

l(i0,][#,](0^|f/(t;,t,)|0,)(0,|p(t,)|0,)(0,|?7(t,,t^)|0;) (19) 



involves two matrix elements of the evolution operator U. To bring out classical features 
more clearly, define 

^ = (</) + <^)/2 (20) 

X = 0-0 (21) 

where (p belongs to the path integral representation of {(f)f\U(tf,ti)\(f)i) and belongs to the 
path integral representation of (0i|?7(tj,t/)|0/). 
Introducing the Wigner functional p^y(v9, vr*^) as 

(y. + x/m^ - x/2) = / M Pwiv, vrl e'^^^ (22) 

we get 

Pw[VN+i,TrN+i] = J[d(po][dxo]pw[^o,T^o] 

N . N . N . N+1 



wwk] I rife] / iiM / nt^' 

k=l k=l •' k=0 •' k= ' 

exp (i 7rfc+i(^fc+i -<^k- TT^+iAt)) 
V fe=o / 



X exp f: Xk (^rr+i - vr^ - E{y,,)At) + 0{x')] (23) 

\ fc=0 / 

Note that without the Wigner transformation of the initial and the final distributions, the 
exponents cannot be arranged in the above finite-difference form that includes the end points. 

The difference between Eq.()23|) and the classical expression Eq.()14j) is the 0{x^) term 
inside the exponential. (There are no terms even in x in this expression because the exponent 
must change sign under the exchange or equivalently, x ~^ ~X-) If we drop the 0{x^) 
terms, then we have exactly the same evolution kernel as the classical theory. However, even 
in this case, this does not mean that pw will be truly classical. If the initial Wigner functional 
contains quantum information, then at any later times will still contain it. 

In perturbation theory, dropping 0{x^) terms means ignoring any Feynman diagrams 
that may contain a loop or loops made up of the off-shell propagators. The validity of such 
operation will be shortly discussed in the next section where we discuss the Feynman rules. 

To summarize this section, the omission of the terms non-linear in x results in the clas- 



sical evolution o: 
and that of Ref. 



the initial distribution. The crucial difference between our analysis here 
is the appearance of the Wigner functional which depends both on the 



coordinate (</?) and the momentum (tt"^). This is important in the perturbation theory 
since the propagators strongly depend on the initial distribution. We therefore discuss the 
propagators and the Feynman rules next. 



III. PROPAGATORS AND FEYNMAN RULES 



Having obtained the classical path integral from a quantum one, we can now derive 
Feynman rules for the classical perturbation theory following the same procedure as in the 
quantum case. To do so, consider the generating functional with the source terms J^(fi — J^x 
and the restrictions that the source terms vanish at the end points ti and tf. Carrying out 
TT^ integrals and xo integral results in 

N 



^1 - ^0 
h 



h 



k=0 



We can then carry out all tt^ integrals to get 



.N+l N 

Z[J^,J^] = / n [^V^fc] / I[[dXk] Pw i^Q, {<Pi - <PQ)/h- E{ipQ)h] 

k=0 k=\ 



X exp 



Xk 

k=l 



h 



N 



E{^k)h + J^h)+iJ2Jk^kh 

k=l ) 



X 



exp X (^-<9V - - ^V'^ ~ "^^j + V 



(24) 



where we switched to the continuum notation for simplicity. 

It is tempting now to develop a perturbation theory from this expression. But we are 
not done yet. Note that the initial distribution depends on the value of the field and its 
time derivative. But the integration is over the field values only. To treat pw[fi, ^i] ^ a 
proper weight we need to transform the measure to include integration over 0j. In classical 
field theory, this can be most easily done by noting that to solve a second order differential 
equation requires specifying cither the values of the field and its first time derivative at the 
initial time or the values of the field at the initial and the final time. Therefore, barring 
caustic points, there is a one-to-one correspondence between the value oi (p at tf and the 
value of (f at ti- The Jacobian of this transformation is 



(25) 
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where S is the action integral. This Jacobian is related to the Van-Vleck determinant 

J[^] = \Det{d^ + m^ + | (26) 

Changing the variable iff to (pi then gives 

Z[J^,J^] = [d7ii][dipi\pw[^Pi,7^i] / 'DipVx 



X J [if] exp 



X 



(27) 



where we have renamed ipi = TCi and there is no integral. One can easily check 

that in the limit Jy, 0, Z ^ 1. This is as it should be since Z reduces to a trace of a 
density operator in this limit. 

Note that without the Jacobian J', the ^ limit will not in general yield 1. This 
is because performing x integral produces 6[{d'^!f + m^Lp + {X/3\)(f^) + J^] which is then 
integrated over (p to yield l/jT". The exact value of the Jacobian, however, depends on how 
one discretizes the equation of motion. This may sounds peculiar, but this phenomenon 
is not new. In stochastic dynamics, it is well known that this kind of functional Jacobian 
depends on Furthermore, the Jacobian can be made 



constant if one choose a particular prescription. This was also noticed in Ref . [ll| . Later in 
this section, we will comment more on the role of the determinant. For now we simply set it 
to a constant and ignore it. As shown shortly, this amounts to the following diagrammatic 
rule: 

(i) Omit diagrams containing tadpoles. 

To get the rest of the Feynman rules, let us examine the two main ingredients of pertur- 
bation theory: Propagators and interaction vertices. The forms of the interaction vertices 
are fully determined by the interaction Lagrangian. In the case of the classical theory, 
the corresponding Feynman rule is: 

(ii) Assign —iX to each Xif^x vertex. 

To determine propagators, one must specify what the un-perturbed state is. In the quantum 
theory, this is the perturbative vacuum which can fluctuate into particle-antiparticle pairs. 
Therefore, the expectation value (O|0(x)0(|/)|O) is non-zero even if (O|0(a;)|O) = 0. On the 
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other hand, in the classical vacuum, both of these quantities are zero because the classical 
vacuum cannot fluctuate. It is literally a state where nothing exists. 

Propagators are determined by the free field limit. Following the derivation in the last 
section, it is easy to see that without the self-interaction, the generating function quickly 
reduces to 

Z[J^,J^\ = j[d'Ki\[dipi\pw[<Pi,T^i\ 

X f VLpjy\5\-d'^Lp-m'^ip + 

regardless of whether it's the quantum case or the classical case. Yet for both cases, the 
combination tp = (0 + 0)/2 satisfies the classical equation of motion whose solution is given 

by 

sni{Ek{t - ti)y 



exp 



i j J^if 



(28) 



+ [ dt'Gretit-t')J^it') 

Jti 



(29) 



Ek I Jti 

The fact that satisfies the classical equation of motion even in the fully quantum case is 
somewhat unexpected. However, it shows explicitly that the form of the initial distribution 
functional plays an essential role in distinguishing the quantum and the classical cases. 

Let's consider the purely classical case first. In equilibrium, the density functional is 
given by the classical Boltzmann factor 



e-^^ci 1 



— exp -/i5:(7r,(^)2 + £;,V,(^)V2 



(30) 



where 7rj(A;) and •piili) are the Fourier components of tt^ and respectively. Using the above 
classical solution and carrying out the and TTj integral then yield 

T /•*/ . /•*/ 



Zc^Jip.J^ = exp 



J2w^ f ' dt f ' dt' J^{t) cos{Ek{t - t'))J^{t') + i I J^Cet 



(31) 



where we have again suppressed the momentum indices for and J^. The resulting Feyn- 
man rules are: 

(iii) In equilibrium, the propagators are 

G7{P) = {T/E,)27:5{p'-m') (32) 
Gthp) = ^1 , (33) 



G^rip) = 
G^^ip) = 



i 

p2 _ _ j^^p^ 



(34) 
(35) 



The same set of rules were also used in Refs. 

mm 

as the classical limit of the 

quantum rules. 

The (fix propagator is the response function which controls generation of classical field 
by a source. Hence it must be there be it classical or quantum {c.f. Eq. (j29|) ). However, in 
the T — > limit, (f(f propagator vanishes due to the fact that classical vacuum does not 
fluctuate. Note also that the classical equilibrium density function is T/Ek just as in the 
case of electromagnetic waves in a cavity. Therefore, if one is to extend this formalism to 
an out-of-equilibrium situation, one should use 

(iii') In non-equilibrium, 

G,7(p, X) = Up, X) 2n 5{p' - m') (36) 

where fc\{p, X) is the density function to be determined and G'^i{p, X) is the Wigner trans- 
formed 2 point function 

Gw{p, X) = j d^r e'^'P^ y) (37) 

with r = X — y and X = {x + y)/2. 

The appearance of the on-shell 5-function in propagator deserves deeper considera- 
tion. Suppose that the initial p[(pi, Xi] is actually independent of Xi- That is, p is diagonal 
in 0. In that case, the Wigner functional contains a 6 functional 5[vrj]. For an illustration, 
suppose 

Pw[^i,^^] = ^exp (^-PY.Elv^{ky/2^ (38) 



In that case, carrying out (pi and vrf integrals results in 

T ( /•* 



Z'{J^, JJ = exp 



(39) 



Since there is no compensating sine term, the ifip Green function no longer exhibits time- 
translation invariance. The ifip correlator is given by 

G'^'^it, t') = ^ cos(Efc(t - U)) cos{Ek{t' - U)) (40) 

Since this is not a function of t — t', this form of the propagator does not conserve energy at 
each vertex. The usual momentum space Feynman diagram techniques will not work. 
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Our ability to use Feynman rules with the usual conservation 4-momentum at each vertex 
depends on whether we get a function only of t — t' for the (fip correlator. In equilibrium, 
the propagators depend only on the difference of the two coordinates. Hence, the total 
4-momentum is conserved at each vertex. In non-equilibrium situations, the propagators 
depend both on the difference r = x — y and the sum X = {x + y)/2. Wigner transformation 
assigns 'momentum' to the difference as done in Eq.()37|). However, unless the X depen- 
dence of Gw{p,X) is slow, the momenta that enter each vertex are not conserved. For the 
derivation of the Boltzmann equation, this slow dependence is therefore essential. 

The above Feynman rules are used to calculate the perturbative corrections in the classical 
field theory. Just as loops appear in the quantum Feynman diagrams, loops appear in the 
classical Feynman diagrams as well ^. This looks like a violation of the theorem that 
states the classical field theory corresponds to the sum of all tree diagrams. However, this 
theorem is derived in the context of the vacuum theory. It does not necessarily apply 
to the in-medium case. Furthermore, the theorem is derived considering only the un-cut 
diagrams where as the in-medium theory necessarily involves cut diagrams (that include (p(p 
propagators). Nevertheless, if this theorem is not to be violated, we must at least show that 
the vacuum loops do not appear in our classical theory. 

That this is indeed the case can be shown as follows. A vacuum loop in a given diagram 
must not contain the density function f^. Since /d appears only in the ip(p propagator, any 
potential vacuum loop diagram must consist entirely of XV' and ipx propagators. Now with 
only the Xxf^ vertex, x is conserved in the sense that it cannot split into multiple x's- Hence 
one can follow x as if it is a fermion. What flows along this line is the time since the ipx 
propagator is a retarded one and the xv^ propagator is an advanced one. Therefore, a loop 
made of (px and/or xyj propagators implies existence of a time- loop which is impossible. It 
follows that there is no vacuum loop contribution in the classical theory if the loop involves 
more than a single propagator. 

The only remaining possibility is a single loop formed by a ipx propagator as shown in 
FiglU which corresponds to 

Sl-loop = ^Gret{t = 0) (41) 

The value of this expression is actually ambiguous as it involves ^(0). This is where the 
Van-Vleck determinant Eq.()26p plays its role. To determine the role of in perturbation 
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X 



FIG. 1: Tadpole diagram in the classical perturbation theory. 

theory, we may exponentiate the determinant using Fadeev-Popov type of ghost fields c and 
c 

Bet{d^ + m^ + J VcVc exp (^i J c{d^ + + {X/2)^^)c^ (42) 

Since there is no in-medium part for the ghost field, the corresponding propagators should 
be 

= -2 W— (43) 

= (44) 

with the rule that a ghost-loop gets an additional factor of (—1). 

First of all, ghost propagators must form a loop since no external line can be a ghost and 
there can be no ghost excitations in the medium. However, since the ghost propagator is a 
retarded one, a ghost loop made of 2 or more ghost lines must vanish. A loop made up of a 
single ghost line is then the only possible non-zero diagram 



Sfrop = -^G',et(t = 0) (45) 



This exactly cancels the tadpole contribution in Eq. ()41|) . Thus the ghost saves the classicality 
of the perturbation theory by removing the ambiguity of defining ^(0). If we define ^(0) = 0, 
then the ghost field contribution disappears and the Van Vleck determinant can be set to 
a constant. Again, as mentioned earlier, this is a known phenomenon. In the study of 
stochastic dynamics, it is well known that how one discretizes can change the form of Van 
Vleck determinant even making it a constant jlfl|. In Ref.^l, since all loops that are formed 
with ifx propagators and/or xf propagators are set to zero, one might say that 6{0) was 
implicitly defined to be zero in that paper. 
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To determine the validity of the above classical perturbation theory as a approximation 
of the underlying quantum theory, let us now consider the quantum case. The form of the 
free field Wigner function can be easily obtained from the ground state wavefunction of the 
simple harmonic oscillator in the Euclidean space (For instance, see Ref. Q].). 

PQb. = ^ exp X: ^^^^^g-^/^^ {n^iky + El^^ikr)^ (46) 

In the small f3 or large T = 1/(3 limit, this goes over to the classical case. Again using the 
classical solution and carrying out the ipi and vTj integral yield 

coth(Efc/?/2) fh rt 



Zq[J'p^Jx\ = exp 



coth^,/j/2j ff^r f ^^^^^^^^ _ ^ ^ j J G,,^j 



(47) 

The propagators are then given by 

G^^ = (1/2 + nBE{Ep))27c6{p^ - m") (48) 
G^^ = (49) 

Gq = (50) 

G^^ = (51) 

where nBE{Ep) = l/{e^p^'^ — 1) is the Bose-Einstein distribution. The corresponding non- 
equilibrium propagator is 

G^J = (l/2 + /Q)2vr5(p2-m^) (52) 

where /q is now the density function to be determined. In Ref. P], these quantum form of 
propagators are used with the classical field Feynman rules (i) and (ii) above. If /g ^ 1/2, 
there is very little difference between this set of quantum propagators and the set of classical 
propagators. However, one of the points made in Ref.j^ was that keeping the 1/2 term in 
Eq.(jn2) makes the classical kinetic equation to match up with the quantum one up to the 
next-to- leading order in /. In view of the fact that this is actually mixing quantum and 
classical descriptions, it deserves a more in-depth study. We will shortly do so in the next 
section. 

Since we have explicit forms of the quantum and the classical propagators now, we can 
discuss the validity of the classical perturbation expansion. The differences between the 
classical and the quantum Feynman rules are as follows (excluding the tadpole diagrams): 
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(i) The quantum i^ip propagator has 1/2 + / whereas the classical one has just /. 

(ii) The quantum interaction includes term whereas this is missing in the classical 
case. 

The first of these items indicates that the classical approximation is valid in the large / 
limit, / ^ 1 so that the appearance of 1/2 does not make a difference. But if / is too big, 




FIG. 2: Ttie first order vertex correction for the classical field. 

then potentially perturbative corrections can be as large or larger than the leading order 
contribution. To be concrete, consider the vertex correction depicted in Fig. |21 Compared 
to the bare vertex, the correction term is smaller by a factor of 0(A/ci). Therefore, the 
classical perturbation theory is valid when /d ^ 1 but also A/d <^ 1. This is the same 
conclusion reached in Ref.jsf but in a heuristic way. 

The next question we should ask is how big the size of the quantum correction is. To 
have an estimate, consider substituting a X(f^x vertex with a Xx^y^ vertex in a classical 
perturbation theory diagram. Since there is no XX propagator, the only way this substitution 
is possible is when the inserted vertices are connected with three (p field and one x field as 
depicted in Fig. 01 In the large / limit, the quantum contribution is 0(1//^) smaller than 
the purely classical contribution. Therefore at a given order in A, the quantum correction 
is always 0(1//^) smaller than the classical one. Again, this is the same conclusion as in 
Ref.Q. 

These estimates also indicate that there is a limit that we can trust the classical per- 
turbation theory depending on the relative size of 1// and A. Suppose that / = 0{l/^/X). 
In the small A limit, this fulfills / ^ 1 as well as A/ ^ 1. However, since 1//^ = 0(A), 
only the leading order classical perturbative correction is reliable. The leading order (in A) 
quantum correction in this case is as big as the second order classical correction. 
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;^(|)^X insertion 




}(^ y} insertion 



FIG. 3: Substituting X^r'x with ^yc"^- 

IV. COMPARISON OF QUANTUM AND CLASSICAL BOLTZMANN EQUA- 
TIONS 

Having derived the Feynman rules, the derivation of Boltzmann equation can now proceed 
as in Ref.j3] for the theory. For convenience, here we show the results from Ref.j5|. The 
collision term of the Boltzmann equation is given by 



2 J (2vr)52^c.(A:0c.(A:,)..(fc3)^^^"^^^ + ^^^^^ " ^^^^^ " ^^^^^^^^^ + - k, - ks) 
i^F{h)F{k2)F{h) + F{p)F{k2)F{k,) - F{p)F{k,)F{k2) - F{p)F{k,)F{ks)^ (53) 

where we used F{p) = /(p) + 1/2 to simplify and combine the terms in Ref.j^f. 

As shown in Section IIIH quantum and classical propagators have different forms. The 
above expression (jS^ uses the quantum form of propagators but otherwise ignores quantum 
corrections. If we are in a purely classical regime, F in the above expression must be changed 
to /. The resulting Boltzmann equation then contains only the terms cubic in /. This is 
appropriate since it is easy to verify that the equilibrium distribution resulting from having 
only the cubic terms is f{E) = T/E. This form leads to the Rayleigh- Jeans catastrophe 
which Planck cured by inventing Quantum Mechanics. 

On the other hand, if we keep the above form, the steady state solution is 



m 



T 
E 



(54) 
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which may be recognized as the first 2 terms of the expansion of the Bose- Einstein distribu- 
tion in the large T hmit. Therefore in the T ^ E hmit, this form is a better approximation of 
the true quantum distribution than the Rayleigh- Jeans form. However unhke the Rayleigh- 
Jeans form, this form of the distribution function becomes negative when T < 2E which 
is clearly unphysical. Therefore mixing the quantum propagators and the classical vertices 
does not necessarily improve the purely classical result. 

Let us now consider why keeping the 1/2 term in the propagator somewhat mysteriously 
reproduce the quantum statistical result up to the next-to-leading order terms in /. Suppose 
we have a scattering process that involves m initial particles (ki) and n final particles (pj). 
The quantum Boltzmann equation must contain the following combination of the density 
functions: 

m n n m 

cq^u fgih) n[i + fQipj)] - n fQiPj) n[i + /qc^^)] (55) 

i=l j=l j=l i=l 

to be consistent with the equilibrium Bose-Einstein density. Expanding in terms of /q's (we 
are in the large /q limit), we have 

n rn n n n m n 

CQ = EiifQiki) n /Qfe)+EE^(^>^)n/Q(^o n /q^) 

s=li=l j=i,j¥=s s=lt=l i=l j=l,jj^s,jyit 

Tfl Tfl TL '11% TfL 71% TL 

-En /Q(^on/Qfe)-EE^(«>^) n /Q(^.)n/Qfe) 

s=l i=l,ijts j=l s=l t=l i=l,i^s,i^t j=l 

+ o(/'"+"-^) (56) 

where we defined ^(condition) = 1 if the condition is met and ^(condition) = otherwise. 
On the other hand, the classical counter-part must contain 

n m n mm n 

^ci = En/ci(^.) n /cife)-E n hmwup,) (s?) 

s=l i=l j=l,j^s s=l i=l,i^s j=l 

to be consistent with the Rayleigh- Jeans form of the equilibrium distribution. It is clear 
that the leading order terms match between the quantum and the classical cases. Now let 

/ci ^ /ci + 1/2 to get 

7% Tfl Tl Tfl Tfl Th 

c^i = E muh) + 1/2] n [/cife) + 1/2] -En luh) + 1/2] n [Up3) + m 

s=li=l j=ljjts s=li=l,ijts j=l 

(58) 

Expanding up to the next order in /d yields 

n m n mm n 

Cci = En/ci(^o n /cife)-E n /cK^on/dfe) 

s=l j=l j=ljjts s=li=l,i^s j=l 
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-inn Tfi n i vn tti tti n 

^ s=l t=l i=l j=lj^s,j=it ^ s=l t=l i=l,i^s,i^t j=l 



+ o(/'"+"-=^) (59) 

which can be recognized to be the same as Eq. ()56|) once we use the fact that 6{s > t) is 
equivalent to 9{s ^ t)/2 under the sums. Thus the prescription, fd — > /ci + 1/2 does produce 
the first two terms in the quantum Boltzmann equation for any general scattering process. 
This also explains why the same kind of consideration also seems to work for hot QCD 17 1 
Although it is curious that this prescription can indeed get the next-to-leading order term in 
/ right, this is irrelevant in getting the right equilibrium limit since that requires retaining 
all the terms in /. 



V. SEPARATION OF HARD AND SOFT MODES 

In view of the discussion in the last section, the better application of the statistical 
mechanics of the classical field is not to get the final equilibrium distribution^ but to use 
it as an intermediate stage where the field degrees of freedom and the particle degrees of 
freedom intermix. In this section, we illustrate how this may be achieved. A more detailed 
analysis of this important problem, however, is out of the scope of the present manuscript 
and will be reported in 

As shown in the last section, the quantum perturbative correction becomes the same size 
as the classical perturbative correction when / ~ 1. As indicated by the equilibrium form 
of /, f{k) becomes small as k becomes larger. Therefore, the classical field description is 
usually appropriate only for the soft modes. 

With the price of introducing a cutoff A (where / becomes ~ 1), the quantum field should 
then be split into the soft and hard modes: 

(p = (Ps + ifh and x = Xs + Xh (60) 

where the soft modes (with the subscript s) have |k| < A and the hard modes (with the 
subscript h) have |k| > A. The soft modes should then be treated as a classical field and 
the hard modes should be kept quantum field. 

^ This was also recognized in Ref.|^ where it was imphed that their formahsm is for the system evolving 
toward equilibrium but not applicable near the kinetic equilibrium. 
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To be concrete, consider the theory with the Hamiltonian given by 
The continuum version of Eq.()23p is 



Pw[<^N+i,'^%+i] = J [di^o][dxo]pw[<^o,'^o] 

X J Vip J Vx J V-n"^ J Vtt^ 



X exp 



(62) 



Separating the hard and the soft modes, the interaction term becomes 



V = ^iXh + Xs)iVh + Vsf + ^iVh + Vs)iXh + Xs? (63) 
Keeping only the terms that contain the most number of ifs yields 

V = ^ (^^Xs^l + ^Xh^'}j (64) 
These terms are linear in Xs and Xh- Hence, the equations of motion are 

aV.-mV.-^y^^O (65) 



and 



d'^cph - m^ifh = ^<^^ (66) 



These equations describe a system where the classical field evolves independent of the par- 
ticles but the quantum mechanical free particles are produces by a classical source. From 
the analogous problem in QED (generation of photons from a classical source), it is an easy 
matter to show that the spectrum of hard modes at the end of evolution is 

ME,M) = ^\JiEkM)\' (67) 

where 

/roc 
d^x dte-*^-''+^^'=*((^,(t,a;))=^ (68) 

This is, of course, in addition to the hard modes which already existed in the initial state. 
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If the cut-off A remains constant at all times, this is the complete solution for the hard 
spectrum provided that we can solve the soft mode classical equation of motion either 
perturbatively or non-perturbatively. In reality, it is not as simple due to the fact that the 
cut-off A should be a functional of ips-, too. In the Color Glass Condensate approach to the 
heavy ion collisions 19], a similar conversion of field degrees of freedom to particle degrees 
of freedom was performed to get partons out of classical non-Abelian gauge field. The line 
of argument given here may provide a firmer ground for such a treatment. One can conceive 
that a kind of 'renormalization group' equation in the manner of Ref. [2^ should exist for 
the cut-off. This will then enable us to describe the whole system in a consistent fashion. 
This topic is currently under investigation. 



VI. SUMMARY AND DISCUSSION 



In this paper, we showed that one can indeed go from the statistical mechanics of quantum 
field theory to the statistical mechanics of particles via classical field theory. Along the way, 
we encountered several subtleties that had to be carefully dealt with. In this regard, we 
improve upon the treatment given in Ref.j^. 

Going from quantum field theory to the classical field theory involves the Wigner func- 
tional which plays the role of initial density function. In other words, the probability of an 
initial state depends on both the value of the field and the value of the first time derivative. 
The form of the initial Wigner functional in turn determines the propagators. In addition, 
our ability to use the Feynman diagrams with energy-momentum conservation at each vertex 
depends very much on the fact that time translation invariance is maintained. For this to 
work at least approximately, the Wigner transformed propagators Gwip^^) (c-/ Eq. ljTTj) ) 
cannot depend strongly on X = (x + ?/)/2. This in turn implies that the density matrix at 
time t must also be approximately of the form 

Inp^El^l + Til^Hu (69) 

where (pk and Tik are Fourier components of the corresponding fields at t and Hk is the 
corresponding free field Hamiltonian. This is a rather strong condition which is in general 
valid only near equilibrium. Whether this can be reconciled with the classical field limit is 
not yet clear and requires further study. 
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Going from the quantum field tlieory to the classical field theory also involves a Jacobian 
in the form of the Van Vleck determinant. The role of this Jacobian turned out to be to 
cancel the single remaining quantum effect still left in the classical perturbation theory in 
the form of the tadpole self-energy. 

It must be also emphasized that to define propagators, one must define what is meant 
by the 'un-perturbed vacuum'. Pure classical vacuum cannot fluctuate as it is literally a 
state where nothing exists. On the other hand, the quantum vacuum fluctuates all the time. 
The importance of this distinction is particularly apparent when one considers the free field 
theory. The form of the path integral for the classical and the quantum theory in this case 
is exactly the same except the form of the initial density functional. To be consistent, the 
classical initial density functional must conform to the classical vacuum property and the 
quantum density functional must conform to the quantum vacuum property. Mixing of the 
two formalism leads to an interesting conclusion that the classical field theory can reproduce 
the results of Quantum statistical Boltzmann equation up to the next leading order in the 
density (Ref.pj). The correct Boltzmann equation obtained from classical perturbation 
theory of course yields the Rayleigh- Jeans form of the equilibrium distribution function as 
it must. 

As for the validity of the classical perturbation theory, in addition to the two conditions 
already mentioned we also argued that there is a limit that a purely classical perturbation 
theory makes sense. For instance, we have shown that if / ~ 1/"\/A, the leading order 
quantum correction is as big as the second order classical correction. And hence, going 
beyond the first order classical perturbation theory doesn't make much sense. Quantum 
effect must be considered after the leading order. 

We would also like comment on the meaning of 'thermalization' here. It is well known that 
an isolated system as a whole cannot thermalize. A simple example is an eigenstate of the 
total Hamiltonian which, by definition, is stationary. This is in contrast to the Boltzmann 
equation where the stationary solution is guaranteed to be the equilibrium distribution. The 
crucial difference of course is whether one is interested in the system as a whole, or just a 
part of the system be it in the momentum space or the coordinate space. It is therefore 
perhaps more natural to apply the current formalism to a particular sector of the system. 
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say the soft modes, and regard the hard spectrum as particle degrees of freedom^. 

In summary, in this paper, we re-examined the derivation of the Boltzmann equation 
from the classical field theory as advocated by Mueller and Son J5|. We pointed out a few 
subtle points in manipulating functional integrals and developing the perturbation theory 
and showed how to deal with them. With our improvement, the framework advocated by 
Mueller and Son has a potential to be a very useful tool in investigating dense many-body 
systems, for instance, the problem of converting field degrees of freedom to particle degrees 
of freedom. Further study in this line is continuing and will be reported in later publications. 
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